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Abstract 

We have previously shown that a ferromagnet-super conductor heterostructure may 
possess a spontaneous current circulation parallel to the interface. This current 
is caused by Andreev bound states in the thin ferromagnetic layer, and can be 
fully spin-polarized. Here we investigate the total energy of the system in cases 
where the current either does or does not flow. We show that the current is a true 
quantum ground state effect, and examine the effect of the current on the different 
contributions to the total energy. 
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The proximity effect between a superconductor and a ferromagnet has been a 
subject of great interest in recent years. Unlike the usual proximity effect be- 
tween a superconductor (S) and a normal metal, a superconductor in contact 
with a ferromagnet (F) is subject to time reversal symmetry breaking which 
splits the up and down spin components of the spin-singlet Cooper pairs. In 
bulk superconducting materials a ferromagnetic exchange splitting leads to 
the FFLO (or LOFF) state in which the superconducting order parameter 
oscillates in space [1,2]. In hybrid S-F structures the same physical effect also 
leads to spatial oscillations [3]. This gives rise to a number of new phenom- 
ena, including: Josephson vr-junction behaviour in S-F-S heterostructures[4,5], 
a giant mutual proximity effect in S-F nanostructures[6], and spin valve [7] 
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Fig. 1. The model system we consider. The ferromagnet (FM) consists of d 
atomic layers {x = —d, ... 0), and is in contact with a semi- infinite superconductor 
{x = 1,2,. . .). A current, J, may flow along y, parallel to the surface. This current 
leads to a magnetic vector potential A in the same direction, and to a B-field in the 
z direction perpendicular to the current, as shown. 

and giant magnetoresistive (GMR) effects [8]. In addition to direct proximity 
effects there are there are also magnetic interaction effects between the B-field 
of the ferromagnet and the Mcissner screening currents of the superconductor. 
These lead to the pinning of superconducting vortices by magnetic domains in 
S-F bilayers[9] and by patterned arrays of ferromagnetic nanodots[10]. More 
complex topologial order parameter textures have also been predicted[ll]. 
Lyuksyutov and Pokrovsky have recently reviewed this rapidly expanding 
literature[12]. 

We have recently examined the S-F proximity effect in the simple model sys- 
tem shown in Fig. 1[13,14,15]. This represents a bulk superconductor with a 
thin ferromagnetic layer on the surface. By solving the Bogoliubov de Gennes 
equations self-consistently for this system we found that the proximity effect 
leads to an oscillatory superconducting pairing amplitude 



sin (x/^f) 



within the ferromagnetic layer. Here Cra is the usual electron annihilation 
operator at lattice site r = {x,y) and C,f = '2t/Eex is the ferromagnetic coher- 
ence length with t the hopping and E^.^ the Stoner exchange splitting within 
the ferromagnet. This oscillatory component has the same origin as the in- 
homogenous FFLO superconducting state, essentially a breaking of time re- 
versal symmetry of the Cooper pairs, and is consistent with experiments such 
as those of Kontos et a/. [3]. The damping of the oscillations as a function 
of distance from the S-F interface is the same as in the usual proximity ef- 
fect with non-magnetic metals, simply the decay of the Cooper pairs in the 
non-superconducting material. 

We also found a much more surprising feature of the system shown in Fig. 1, 
namely a spontaneous spin-polarized current, J, parallel to the S-F interface[13,14,15]. 
An unexpected feature of this current is that it switches on or off suddenly as 
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a function of the number of ferromagnetic layers, d. divided by the coherence 
length ^p. In fact we found the the stable self-consistent solution of the Bogoli- 
ubov de Gennes equations has a finite current precisely when the dimensionless 
parameter, © = 2.79dEgx/ {nt) is equal to an odd integer, n = 1,3,5.... This 
surprising condition precisely corresponds to the conditions for the supercon- 
ducting order parameter at the top ferromagnetic layer rf, 0) to be zero. 
As show in [15] the top layer order parameter ci.O) oscillates as a function 
of Egx, changing sign whenever is an odd integer. Further analysis showed 
that this unusual condition for the current carrying solution arises from the 
spectrum of the Andreev bound states within the ferromagnetic layer. An in- 
stability arises whenever a pair of Andreev bound states become degenerate 
at the Fermi level. The self-consistent solution with the current flow splits the 
degenerate pair, and leads to a spontaneous parallel current flow, J along the 
S-F interface. 

In this paper we investigate in more depth this surprising sudden turning on 
and off of the current carrying solution. We calculate the total energy of the 
system shown in Fig. 1 both with and without the current flow. This allows 
us to explain more clearly the physical origin of the current, and to confirm 
that the current is indeed a true quantum ground state effect. In the next 
sections we first outline the formalism and computational methodology and 
for our calculations, before we then present the results. 

Our model Hamiltonian is given by the negative U Hubbard model. 
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where, the nearest neighbour hopping integrals are tij = —te ''^^''^ in 
the presence of a vector potential A{r). The site energies £i(j are on the 
superconducting side and equal to ^E^^a on the ferromagnetic side, fi is the 
chemical potential, and C/j is C/5 < in the superconductor and zero elsewhere. 
Here c^, (cia) are the usual electron operators and fiia — cf^Cia- 

We solve the above model in the Spin-Polarized-Hartrec-Fock-Gorkov {SP- 
HFG) approximation. We work in the Landau gauge where B = (0, 0, Bz{x)) 

— * 

and hence A — {0, Ay{x),0) (see Fig.l). Furthermore, we assume that the 
effective SPHFG Hamiltonian is periodic in the direction parallel to the in- 

— * 

terface and therefore we work in k space in the y direction but in real space 
in the x-direction. As usual, self-consistency is assured by the relation: 
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Fig. 2. The difference in total energy between solutions with and without current 
flow. The difference is large whenever the parameter = 2.79dEex/{'T^t) is an odd 
integer, and the current flow is stable. For other values of the energy difference 
is zero, and no current flow occurs. 

where G°!^{uj,ky) is the 4x4 Nambu Green function and f{uj) the Fermi 
function. 

Since our model includes a vector potential the solution may imply a non-zero 
current. For spin up electrons, in the ?/-direction this can be calculated from: 

J,ta)(n) = -2etJ2s^niky - eA,(n)) / duj-lmGlf'\uj, ky)f{uj). (4) 

h J TV 



The current will give rise to a vector potential A{f) which is obtained solving 
Ampere's law, ^ ^^i^'' = — 47rJy(x). This full set of equations can be solved self- 
consistently, using the method described in [15]. The various contributions to 
the total free energy can also be calculated, and the expressions for the various 
contributing terms are given in [16]. 

The coupled set of non-linear equations given above can be solved self-consistently 
both with and without a finite current flow, J. Fig. 2 shows the change in to- 
tal energy of the system in Fig. 1 between the fully optimized self-consistent 
solution and the one where the current is constrained to be zero. We used the 
parameter values: d = 20 layers, temperature T = 0.01 and U = —2.345 in 
units of the hoping integral t. Comparing the total energy of these two cases 
we can see, in Fig. 2 that the current carrying solution is indeed the one of 
lower total energy, when the current carrying state is found. These are the 
series of sharp minima, each one occurring where the dimensionless parameter 
B (defined above) is an odd integer, and where the pairing order parameter 
x(r) changes sign on the topmost layer of the ferromagnet. 

In between these sharp minima, the energy difference goes to zero, because 
when 9 is away from odd integer values the ground state has no spontaneous 
current flow. One can also see in Fig. 2 that the magnitude of the energy 
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Fig. 3. Two contributions to the total energy difference, (a) Left, shows the change 
in electron kinetic energy Ekin- (b) Right shows the change in energy from the 
current term Ej. 



difference gradually declines as the exchange splitting Eex becomes larger. 
This is because the oscillatory x(r) in the ferromagnet (Eq. 1) becomes more 
and more heavily damped as the exchange splitting increases. The energy gain 
associated with the current declines, and for very large E^x the spontaneous 
current J is suppressed. The modulations in the size of the energy gain shown 
in Fig. 2 are due to oscillations in the amplitudes of the Andreev states in the 
density of states at the Fermi energy. 

One can analyze these contributions further by examining the various terms 
contributing to the total energy. In Fig. 3 we show the change in electron 
kinetic energy (tijcf^Cjo)- Clearly this is very similar in magnitude and form 
to the total energy change of Fig. 2, showing that the spontaneous current 
is mainly driven by the desire to lower the electron kinetic energy. Other 
contributions to the total energy change, such as the 'gap' term, /Ui, 
shown on the right in Fig. 3 are not only smaller in magnitude, but can also 
be either positive or negative. This shows that the spontaneous current flow 
has only a marginal effect on these contributions to the total energy. 

In conclusion we note that these results are consistent with the picture devel- 
oped in Ref. [15]. The spontaneous current flow is generated to lift a Fermi 
surface degeneracy between Andreev bound states, and is primarily driven by 
kinetic energy optimization among the Andreev bound states. The conditions 
for this current to be observable experimentally are (i) that the ferromagnetic 
exchange E^.^ be not too large (i.e. a weak ferromagnet), (ii) the ferromag- 
netic thin fllm thickness should be thinner than the electron mean free path, 
d < I [17], and (iii) the F-S interface transparency must be greater than 
about T] > 0.3 [17]. In systems satisfying these conditions, the existence of 
the spontaneous current may be conflrmed experimentally by observing the 
related magnetic flux ~ 0.0125$o per plaquette, or by STM spectroscopy on 
the Andreev states near to the Fermi level[15]. 
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